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Abstract 
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1 Introduction 



Cold ions confined in electromagnètic traps and cooled by means of làsers 
are a very important experimental system both for the study of fundamen- 
tal physics, such as cold non-neutral plasmas or quantum dynamics, and for 
technological applications such as optical frequency standards. In the past 
few years systems of this kind have been the subject of intense study as a 
possible architecture for the realization of a quantum computer [1, 2, 3] (for 
reviews of progress towards this goal, see [4, 5, 6, 7]). In comparison with 
other experimental systems of investigations of fundamental quantum phe- 
nomena, cold-trapped ions offer a number of advantages. In particular, it is 
a relatively "clean" system whose behavior is well characterized by theory; 
simplified models of interactions between two-level systems and quantum har- 
mònic oscillators can be realized experiment ally Considerable achievements 
in the production of non-classical states of matter [8] and other fundamental 
physics problems such as high-efficiency measurements of Bell's inequalities 
[9] or tests of cavity quantum-electrodynamics [10] have been made. 

As is well known [11], it is impossible to confine charged partides by 
electrostatic forces alone. To overcome this problem, the radio-frequency 
Paul trap was developed: such devices create an effective binding potential 
in the x and y directions while a weaker static potential is applied in the 
z direction [12]. When two or more ions are confined in such a trap, they 
will repèl each other due to the Coulomb force, resulting in confined charged 
plasmas with very low densities. When sufficiently cold, the plasma will 
condense into a crystalline state. In the highly anisotropic traps used for 
some atòmic clocks [13] and for quantum computing, this crystalline state is, 
for small enough numbers of ions, a simple chain of ions lying in a straight 
line along the axis of weak binding. In what follows, we shall be examining 
the oscillatory of ions confined in such a linear configuration. 

The ions are assumed to be sufficiently cold that they undergo only small 
oscillations around their equilibrium positions. In this case, the Lagrangian 
of the motion can be expanded as a Taylor series around these equilibrium 
positions. As we shall see, the first term of this expansion is a constant 
(which has no effect on the motion), the second term vanishes (by virtue 
of the condition for equilibrium); the third term describes coupling of one 
ion's motion to all of the others (and can be easily resolved into a series 
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of normal modes describing the ions' collective oscillations); the fourth term 
(whose consequences are the main theme of this paper), describes an intrinsic 
coupling between the different modes. 

The effects of non- linear ities in the coupling of harmonically bound par- 
tides is a problem with a long history: for example Fermi, Pasta and Ulam 
investigated such dynamics computationally in 1954-55, using the early Los 
Alamos MANÍAC computer, with a view to explaining the equipartition of 
energy between modes as the system reaches thermal equilibrium [14]. In 
the quantum realm, the non-linearities lead to couplings between the nor- 
mal modes of oscillation. In certain ways this coupling is analogous to 
non- linear optical effects [15], in that, provided certain conditions are met, 
quanta of one oscillatory mode can be down-converted into twinned pairs of 
quanta of other modes. Obviously, since the Taylor expansion is an infinite 
series, higher terms, describing more complex multi-mode couplings, will be 
present; we will, however, not consider them in this current work. 

The paper is organized as follows. In Section 2 the Lagrangian of the 
ion motion is presented and the derivation of the normal modes of the ion 
oscillations is reprised; the classical description of mode-mode coupling in 
terms of the Hamiltonian, is formulated. In section 3 the properties of the 
mode-mode coupling coefhcients are discussed, and section 4 discusses the 
quantum theory of this coupling and the population transfers that it can 
cause. Section 5 summarizes our conclusions. 

2 Lagrangian of the ion motion in the trap 

Consider N identical ions of mass M and charge Q confined in an effective 
three-dimensional harmònic potential 1 . The position of the n-th ion will 
be denoted (x n i,x n 2,x n3 ) (Fig.l), where the ions are numbered in order of 
increasing value of their axial positions x n3 , so that n > m implies that 
Xn3 > x m3 . Besides the trap potential, each ion experiences a Coulomb 
interaction with each of the other ions. Thus the total potential energy of 

l· In reality, the ions are experiencing a timc varying force whose effects can be mod- 
eled as an effective harmònic motion plus a high-frequeney oscillation called micromotion. 
Experimental techniques for minimizing this effect are explained in ref.[16]. Here we will 
assume that the amplitudes of transverse oscillations are sufficicntly small that micromo- 
tion can be neglected in this paper. 
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the ions is given by the following expression: 

N 3 



2 ±^é^ 87re 



ra=l j=l 



N 

E 

n,m=l 



3 
i=l 



(2.1) 



where u;i, c^2 and are the angular freqüències of the trap in the three 
directions. For simplicity, we will assume that the trapping potentials are 
equal in two transverse directions (1 and 2) and that the trapping potential 
in the axial direction (3) is much weaker i.e. 



(2.2) 



where a < 1 is a dimensionless parameter characterizing the anisotropy of 
the trap. 



x — X 1 




equilibri um position of n-th ion 
/ (0-0- f«„) 



position vector 
of n-th ion, (x nP x n2 , x n3 ) 



axial direction, z =x 3 

displacement 

of n-th ion, (^, ^ Ç n3 ) 



Figure 1: Illustrating the notation used in this paper. 

The equilibrium positions of the ions, denoted by x n i, (n = 1, . . . , N), 
(i = 1,2,3) are determined by the following equations: 



dV 



= (n = l,...,N) (i = 1,2, 3), 



(2.3) 



where the subscript denotes that the partial derivatives are evaluated 
at x ni = x n i. The solutions of these equilibrium equations in strongly 
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anisotropic trapping conditions (so that the ions are aligned along the weak 
axial direction 3, i.e. x ní = x n2 = 0) have been investigated by various 
authors [17, 18, 19] (see also [20, 21] for analysis of ion crystals with two 
different species of ion). It will be convenient to use the dimensionless equi- 
librium positions defined by u n = x n s/£ where the length scale í is defined 
by 



1/3 



\ 47re Mcj| 



(2.4) 



The displacements of the ions from their equilibrium position are de- 
noted £ ni , i.e. 

The Lagrangian describing the motion is then 



L = T-V 

M N 3 
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Making a Taylor expansion about the equilibrium positions, this may be 
approximated as 
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(2.7) 



Neglecting both the constant V (which has no effect on the dynamics) and 
the higher order terms and evaluating the partial derivatives explicitly (see 
Appendix A), the Lagrangian may be approximated by following expression 
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The first term represents the ions oscillations along the axial (i = 3) direction, 
the second oscillations in the transverse directions (i = 1, 2) and the third 
term represents coupling between these oscillations which are a direct and 
unavoidable consequence of the Coulomb interaction between the ions. The 
tensors A mn , B mn and C mnp are given by 



A, 



TV 
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p=ém 
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if m = n , 



if m 7^ n, 



(2.9) 



(2.10) 
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(2.11) 



where sgn(x) stands for the sign of x and 5 mn is the Kronecker delta. The 
third rank tensor C mnp is symmetric under any exchange of two subscripts; 
thus, for example, C 2 24 = C* 2 42 = C 42 2 = — l/(ti2 — M 4) 4 while C442 = C244 = 
C 42 4 = l/(«2 — w 4 ) 4 . All the elements of A mn and C mnp can be calculated 
numerically from the dimensionless equilibrium positions of the ions. 

The matrix A mn is real, symmetric and positive definite. Thus its eigen- 
values fi p are non-negative. The eigenvectors bffl (p — 1, . . . , N) are defined 
by the following formula: 
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(2.12) 



where the eigenvectors are numbered in order of increasing eigenvalue. The 
eigenvectors form a complete basis so that 
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(2.13) 
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For any value of N, the two first eigenvectors are the center of mass mode 

= l/y/N(l,l,...,l), (pi = 1 ) and the stretch mode b {2) = (l/jV)(wi, « 2 , . . . , u N ) 
(fx 2 = 3) where M = \Jj2n=i u n- AU the eigenvalues and eigenvectors can 
be calculated numerically; their approximate vàlues are given for example 
in [18] for 2 to 10 ions. The dynamics of these modes have been confirmed 
directly in experiments [22]. From eq.(2.10) we see that B mn has the same 
eigenvectors as A mn but different eigenvalues: 



N 



]T B mn b® = lp bt ] (p = l,...,iV), (2.14) 

n=l 

where the eigenvalue 7 P is related to the longitudinal eigenvalue by the for- 
mula 

7 P = - + - " — • (2-15) 
lp a 2 2 v ; 

If a > a crit = 2/(//jv — 1) the matrix £> nm is not positive defmite; this im- 

plies unstable transverse oscillation modes, and results in the formation of 

"zig-zag" crystal structures [23]. We shall assume however that ct is kept 

sufficiently small that this situation does not arise. Note that, with this con- 

vention, the transverse oscillation modes are numbered in order of decreasing 

eigenvalue, so that the p — 1 center of mass mode has the largest eigenvalue. 

Using these eigenvectors, the normal modes of the ions' oscillations in 

the three directions are defined as follows: 

N 



*p(t) = £^nl(í) , 
n=l 
N 

Y P (t) = Y.WUit) , (2.16) 

n=l 
N 

z P (t) = E Wu(t) . 
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The Lagrangian eq.(2.8) can be rewritten in terms of these normal modes 
quite straightforwardly. Not surprisingly, one finds that, neglecting the 
fourth order term, the Lagrangian becomes a sum of Lagrangians of harmònic 
oscillators. The canonical momenta conjugate to X n is Tl* = dL/dX n (with 
analogous definitions for IT^ and n*). Using these momenta, the Hamiltonian 
for the ion motion is 

H = H + H I: (2.17) 
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where H is the Hamiltonian for all of the uncoupled collective oscillations, 
i.e. 



H = 
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(2.18) 



z/ n = ^y^ün being the angular frequency of the n-th mode in the z direction 
and íl n = the angular frequency of the n-th transverse mode. The 

Hamiltonian Hi describes the perturbation which, in certain circumstances, 
couples these various modes. It is given by 

Meu 2 N 

Hj = z ^2 D mnp Z p (2Z m Z n — 3X m X n — 3Y m Y n ), (2-19) 

m,n,p=l 

where the mode-mode coupling coefficients D pqr are defined by 

N 

D pqr = E C^b^. (2.20) 

l,m,n=l 

Note that the sign of D pqr depends on the sign of the eigenvector h^ p \ which 
is undetermined. We will adopt the convention that the sign of the iV-th 
component of b^ p ^ is always positive, for all modes and vàlues of N. 

We should mention that other mechanisms by which resonant interac- 
tions between ion motion in different directions are coupled have been stud- 
ied. In particular Werth et al. at the University of Mainz have investigated 
experimentally the parametric coupling of different stable configurations of 
molecular ion clouds [24]. Cross-mode couplings due to static field imperfec- 
tions, which can be a cause of heating or of decoherence, were also discussed 
in ref.[4] sec 4.1.8. 



3 Properties of the Coupling Coefficients 

In this section we discuss some of the properties of the coefficients C mnp and 
D m np, which are central to the subsequent development of phonon-phonon 



8 



interactions in ion traps. The first important symmetry property stems di- 
rectly from the definition of C mnp , eq.(2.11): The tensor C mnp is symmetric 
under exchange of two subscripts; further, because of the definition (2.20), 
so is D mnp . From the definition of C mnp eq.(2.11), we also have the property 
Cmmn = —Cnnm (provided that m ^ n) . These properties imply that 

N 

E Crnnp = 0. (3-1) 
p=l 

This can be shown as follows: if m ^ n, we have Y, p =i C mnp = Cmnn + 
Cmnm, since C mnp is always zero if all three indices are different. However, 
since C mnm = C mmn = -C nnm = -C mnn , the sum is zero. If m = n, 
T,p=iC mmp = C mmm + T,p=i, P jÈ m C mmp , which, by the definition of C mmm , 
eq.(2.11), must also be zero. The identity eq.(3.1) implies that, since 

(i/Vïv){i,i,...i}, 

i TV N 

D mnl = -±= E t^f E C m = 0. (3.2) 

Thus, when considering the three-mode mixing described by eq.(2.19) the 
center of mass (p = 1) mode has zero coupling strength to any other mode. 
Physically, this can be explained as follows. The p — 1 mode is special since 
it is the only mode in which the center of mass of the crystal is displaced: 
as can be seen from eq.(2.13), J2m=i brn = if p 7^ 1. Thus to excite or de- 
excite the p — 1 mode requires the application of an external force to change 
the momentum of the crystal as a whole, while no such force is required to 
transfer energy between any of the other modes. 

This result has important consequences: firstly, if one wanted to avoid 
mode-mode coupling entirely, for example when performing quantum lògic 
operations between ions, it would behoove one to use the center of mass 
mode as the "quantum bus" as originally proposed by Cirac and Zoller [1]. 
However, in some experiments [25], it is found that the quantum state of 
the center of mass mode is rapidly degraded by excitations from various 
extraneous sources. Our result eq.(3.2) tends to rule out coupling to other ion 
oscillation modes which may have been imperfectly cooled as a sources of such 
heating. Furthermore, when the center of mass mode is heated preferent ially, 
it will not cause heating of the other modes. 
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A second interesting property of C mnp can be expressed as follows: 

N 1 

y ' Up C mnp — [S mn A mn \ , (3-3) 
P =i 1 

where A mn is the ion coupling tensor defined by eq.(2.9). This result can be 



mnmU'm 



demonstrated as follows: if m ^ n, then X) p =i C mnp u p = C mnn u n +C, 
{u n — u m )C mnn = l/\u n — u m \ 3 , where the last step results from the definition 
of C mnn in terms of the scaled equilibrium positions u p . For the case m = n, 

u m \ 3 . Comparing these expressions involving the -u p 's with the definition of 
A mn , eq.(2.9), we obtain eq.(3.3). This result implies the following 

N N Uk 1 — /i 

D r nn2 = ^ V T? = OA/^" 1 ™' 

i,j=l k=l JV ZyV 



where, as before, J\í = yY^n=i u n- Thus the coupling of the first stretch mode 
(p—2) is constrained so that it will only be coupled to a single other mode 
rather than to two distinct modes. 

For N = 2 and N = 3, we can determine the coefficients D mnp alge- 
braically, using the exact expressions for the equilibrium positions and mode 
vectors which can be obtained in those two simple cases. We find that 

N = 2:D 222 = -2s, 

n^ = -±(t^ 



N = 3:D 233 = --^y .^ = -^y • (3.5) 

All the other coefficients must be determined numerically, although this is a 
reasonably straightforward task. The approximate numerical vàlues of the 
non-zero coefficients for 2 to 10 ions are given in Tables 2a and 2b (Appendix 
B). 

4 Quantum motion of the ions 
4.1 Quantization of the Hamiltonian 

We can now consider the quantum motion of the ions by introducing the 
following position and momentum operators 



z n (t) - z B = y ^(an + ai) K(t) - K = i)J^^(àl - à n ), 
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x n (t) — > X n 

Y n (t) - Y n 
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-(k + bi) K(t)^fi x n 
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'Mh a 



(bí-k), 



(4.1) 



where, as before, the Z-coordinate refers to motion along the axial direction 
of the trap, while the X- and F-coordinates correspond to displacements in 
the transverse direction. These operators obey the canonical commutation 
relations 



x p , 

The Hamiltonian becomes: 
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H = H + Hj, 
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where 
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(4.4) 



where e is a dimensionless quantity which characterizes the strength of the 
non-linearity; it is given by 

n i/6 

(4.5) 



£ = 



4v/2 



^ 3 



L«fsc 



Mc 2 



where af sc is the fine structure constant. Alternatively, e = a/4£, where 
a = ^Jh/2Müü 3 is the approximate width of the wavepackets of the individual 
ions, and £, defined by eq.(2.4), is the length scale of the order of the ions' 
separation. divided by the ions' separation. Since both D mnp and fi p are, at 
most, of the order of unity, while 7 P is assumed to be much larger than unity, 
we have: 



H 



e. 



(4.6) 
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Vàlues of e for a variety of ions and traps are given in Table 1. These show 
this parameter tends to have a small value in experimental circumstances, 
implying that the treatment of Hi as a perturbation can be justified. 



Ion 


u 3 


€ 


EU 3 


9 Be+ 
40 Ca+ 

88g r + 

112 Cd + 


(2?r) 5.0 MHz 
(2?r) 2.0 MHz 
(2?r) 200 kHz 
(2vr) 2.8 MHz 


1.06 xl0~ 3 
7.09 xl0~ 4 
4.24 xl0~ 4 
6.32 xl0~ 4 


(2tt) 5.30 kHz 
(2tt) 1.42 kHz 
(2?r) 85 Hz 
(2vr) 1.77 kHz 



Table 1: Vàlues of the parameter e for various ions and trapping freqüències. 
We will denote the eigenstates and eigenvalues of H as: 

|tf> = \{n* 1 ,...,n%},{nl...,n y N },{n z 1 ,...,n z N }) (4.7) 

N 

Eia = hJ2n p (n; + n y p )+is p n z p - (4-8) 
P =i 

4.2 Resonance conditions 

Fermi's golden rule [26] implies that population transfer between two eigen- 
states \i) and |/) of Hç, induced by the interaction Hj can only occur if the 
matrix element (f\Hj\i) is non-zero and if the energies of \i) and |/) are equal 
2 . These conditions allow us to neglect a considerable number of the terms 
that occur when eq.(4.4) is expanded. For example, because the square roots 
of the eigenvalues of A m rí are all irrational numbers (with the sole exception 
of /ii), it follows that ± (y/^i yfjï^± v^ün) 7^ 0; hence we can neglect all 
the terms of the form (à p + àp) (à m + à^) (à n + a^). Similarly terms involving 
three creation or three annihilation operators can be ignored. Making these 
approximations, and using the symmetry properties of D mnp and the com- 
mutation relations (4.2), we obtain the following simplified expression for the 
interaction Hamiltonian, (in the interaction picture) 



N 

Hj tu —3ekuj, x " ' 



3 rr- — — — - 

m,n,p=l \f~1rrTinHp 



D, ,,,,,, . /;-:-; n iA ( m ï p üj z t 



2à p (bl.bn + c^c n ) e l 



2 This is equivalent to the energy matching condition in parametric down-conversion. 
Since here the modes do not represent travelling waves, there is no momentum conservation 
(or phase matching) condition as there is the non-linear òptics. 
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+ h.a., 



(4.9) 



where h.a. stands for the Hermitian adjoint of the preceding term, and 
^■mnp = sfïm ± \^fn — ^[Wp- We shall refer to the term proportional to 
exp(iA^ p u z t), which involves the creation of a transverse phonon and si- 
multaneous annihilation of both a longitudinal and a transverse phonon (and 
the reverse process, contained within the hermitin adjoint part of eq.(4.9)), 
as a resonance of the first kind; the term with exp(i A^ p u; z t) , (which cre- 
ates two transverse phonons in different modes by annihilating a longitudinal 
phonon) will be called a resonance of the second kind. 

Since the only experiment ally controllable parameter in the definitions 
of is the anisotropy parameter a, it is natural to ask at what vàlues 

of a a resonance can occur. A necessary condition for resonance (i.e. for 

0) is as follows: 



A(+) 

mnp 



or A^ 



a 



16 n P 



(4.10) 



Using the known vàlues of the eigenavalues \i v (ref. [18], table 2), appropri- 
ate vàlues for a can be found straightforwardly. However, not all of these 
vàlues of a given by (4.10) correspond to a resonance, and one must deter- 
mine whether or not the appropriate condition for either type of resonance 
is satisfied by direct substitution into the formulas for A^ p - Since ^fJTp ^ 0, 
resonances of the first kind are only possible when two distinct transverse 
modes are involved (i.e. m ^ n); resonances of the second kind can occur 
involving a single transverse mode (i.e. m can be equal to n). Resonances of 
the first kind tend to be very weak; the coupling coefficients do not exceed 
1.6 x 10~ 3 for L < 10, and are always zero for L < 6. The vàlues of a for 
the two types of resonance are given in Tables 2a and 2b (Appendix B) and 
are plotted for different numbers of ions in figure 2. Only those vàlues less 
than the critical value a cr ü are included, since our analysis is based on the 
assumption that the ions are in a linear configuration. 

One important result stems from this analysis. Resonance cannot occur 
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Figure 2: Vàlues of the anisotropy parameter a for which resonant three- 
mode mixing can occur. The horizontal/verticle crosses represent resonances 
of the first kind (A^ = 0) and the diagonal crosses represent resonances 
of the second kind (AW = 0). The upper limiting curve represents the 
critical value a crit , above which the linear confirguration of the ions becomes 
unstable. The lower curve represents the highest value of a for which resonant 
mode-mode coupling cannot occur (see eq.(4.11)). 



for vàlues of a below the minimum value given by 

16/íjv 4 



a 7 



Un (9/xjv + 2/íat-i - 8) + 
4 

2>ijl n — 2 



> 



3/ijv — 2 



if N is odd, 



if is even. 



(4.11) 

Thus, should one wish to avoid any form of mode-mode coupling due to the 
forth order terms we are considering in this paper, this can be ensured by 
using a trap with an anisotropy parameter smaller than 4/(3/íjv ~ 2). This 
can also be expressed in terms of the critical anisoptry for onset of zig-zag, 
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a cr i t : there will be no three-mode mixing if a < 4a cr i t / (a cr it + 6). 



4.3 Example of mode-mode coupling 

As a specific example of the above general analysis, let us assume that we 
have six ions confined in a trap with anisotropy factor a—O. 09151. For this 
value of a, a resonance of the second kind (A^j = 0) occurs for the three 
modes (m,n,p) = (5,6,5) and also for (m,n,p) = (6,5,5) (since there is 
symmetry between m and n for resonances of the second kind). 

The coupling Hamiltonian can be simplified by neglecting all off-resonant 
terms; for this case we get 

~ -3e?iu} 3 — ^^=à 5 (blbl + c£c 6 ) 

^57576 V 1 

+ - 3ehcü 3 — - 655 q 5 (UU + 44) + h.a. 

^Ü57675 v > 

= -ehüü 3 a 6j ° 565 q 5 (blbl + clcl) + h.a., (4.12) 

^57576 V ' 

where we have used the fact that D 655 = D 565 and that 6 5 ,6g =0 (and 
similarly for the cjj. From Table 2 we see that -D 655 =4.2528 and a=0. 09151; 
from Table 2 of ref.[18] we find /i 5 = 13.51; using eq.(2.15) we find 75 = 4.6709 
and 7 6 =2.2949; thus 6L> 655 / ^57576 = 7.3556. 

As an example we will consider the case of coupling between three states 
defined by 

|*> = | {0,0, 0,0, 0,0}, {0,0, 0,0, 0,0}, {0,0, 0,0, 1,0}), 

|$) = I {0,0, 0,0, 0,0}, {0,0, 0,0, 1,1}, {0,0, 0,0, 0,0}), (4.13) 

Ix) = I {0,0, 0,0, 1,1}, {0,0, 0,0, 0,0}, {0,0, 0,0, 0,0}), 

where we have used the notation for the motional quantum states of the ion 
modes defined in eq.(4.7). Note that the coupling defined by eq.(4.12) only 
involves the 5th longitudinal mode and the 5th and 6th transverse modes; 
the states of the other modes are unaffected; thus although we have set these 
other modes to have zero population in the definitions of \<p) and \x), 
their state can be arbitrary without changing the analysis. The interaction 
picture state describing the system is given by: 

i*(f)> = mm + «KOI*) + x(t)\x)- (4.i4) 
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Since, due to the resonance condition, H , Hj — and the interaction 
picture coupling Hamiltonian is: 

eï Èot H ie -ï Èot = Hj. (4.15) 

Thus the Schròdinger equation describing the evolution of the system is given 
by 

ihl\* I (t)) = H I \* I (t)), (4.16) 

which implies that the probability amplitudes obey the following set of cou- 
pled equations: 

#(í) = -r[0(í) + x (t)], 

i<j>(t) = -r^(í), (4.17) 

ix(t) = -r^(í), 

where T = ecj 3 (6-D 565 / • The solution of these equations can be 
obtained straightforwardly using Laplace transforms: 

m = m cos (Vlrt) + i m ~t x(0)) sin (Vlrt) , 

0(í) = ^5 sin (v^H) + 0(0) cos 2 (IQ - X (0) sin 2 ^ ,(4.18) 

X(í) = ^sin(V2rt)-0(O)sin 2 ^ +X (O)cos 2 ^. 

These solutions are plotted in Figure 3 for the case ip(0) = 1 (i.e. the 
longitudinal phonon mode initially excited) and in Figure 4 for the case 
0(0) = 1 (i.e. two transverse modes in the x-direction initially excited). 
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TIME (scaled units) 

Figure 3: The evolution of the populations when the initial state is \ip). The 
thin line represents |t/>(£)| 2 and the wide line |0(í)| 2 = |xW| 2 - Time is in 
units of I/T. 

4.4 Entanglement 

Let us consider in detail the first example, shown in Fig.3. The crystal of six 
ions was prepared with a single quanta in the fifth oscillatory mode, i.e. the 
state 

| (0)> = | {0, 0, 0, 0, 0, 0} , {0, 0, 0, 0, 0, 0} , {0, 0, 0, 0, 1, 0}). (4.19) 
After a time t\ = n/2^/2T the crystal will have evolved into the state 

Mh)) = -^(| {o,o, o,o, o,o}, {o,o, o,o, i,i}, {o,o, o,o, o,o}) 

+ | {0, 0, 0, 0, 1, 1} , {0, 0, 0, 0, 0, 0} , {0, 0, 0, 0, 0, 0})) . (4.20) 

This later state has the property of being entangled: it cannot be written as 
a product of the state of 5-th and of the 6-th oscillatory modes. This can be 
thought of as analogous to the type of polarization entangled states of pho- 
ton pairs generated in optical down-conversion experiments [27]. The two 
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TIME (scaled units) 

Figure 4: The evolution of the populations when the initial state is \<f>). The 
thin line represents |^(í)| 2 , the dashed line |0(t)| 2 and the wide line |x(í)| 2 . 
Time is in units of 1/T. 

transverse directions of oscillation are analogous to the photon polarizations, 
while the oscillatory mode is analogous to the spatial mode of the photons. 
The initially excited longitudinal oscillations play the role of the pump làser. 
Controlling this process in ion traps will be achieved by switching the trap 
anisotropy parameter a into the appropriate resonance parameter for the ap- 
propriate amount of time, then rapidly switching it to a non-resonant value: 
this will have the effect of turning on the phonon-phonon interaction for a 
specific time. Furthermore, the technology of preparing oscillatory quantum 
states of trapped ions is at the moment more advanced than current state 
of the art for preparing quantum states of the electromagnètic field. Specif- 
ically, one can hope to prepare the "pump" longitudinal mode in a single 
excitation Fock state, while creating single photons is still somewhat prob- 
lemàtic. Thus in principle, deterministic preparation of oscillatory entangled 
states may be attained, in contrast to the non-deterministic, post-selection 
state preparation in quantum òptics experiments. Furthermore the method 
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of entangled oscillatory state preparation described here, combined with the 
variety of techniques for creating entanglement of the infernal degrees of 
freedom raises the possibility of creating hyper- entangled states in ion traps 



In this paper we have analyzed the dynamics of ions confined in harmònic 
traps. The normal modes of the ions' oscillations are intrinsically coupled 
by the Coulomb interaction. This "three-mode mixing" is in some ways 
analogous to non-linear optical effects. There are two important results which 
our analysis has revealed: first, if one wishes to avoid these effects altogether, 
the trap anisotropy, characterized by the parameter a defined by eq.(2.2), 
must be less than a certain value given by eq.(4.11). Secondly, this effect in 
principle may be exploited to create entangled motional states of the ions, 
by carefully controlling the ions' dynamics, as outline in section 4.3 and 4.4. 
How the resulting entangled states may be exploited for tasks in quantum 
information is a promising avenue for further investigation 
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Appendix A: Derivation of Eqs.(2.9), (2.10) and 



In this Appendix, we show how to derive the Lagrangian eq.(2.8) from the 
eq.(2.7). We will use the following notation: 



[28]. 



5 Conclusion 




(mn) 



(A.l) 



nií 
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i=l 



(A.2) 



The potential energy of the ions is 



M 



N 3 



2^2 N 



Z 2 e 



^ = tEE« + ^- E 



n=l i=l 



87re ^ , i?( mn ^ 



(A.3) 



Thus we have: 



<9x r 



e 3 uj " Rt p) 



(A.4) 



<9 2 \/ M ^ 2 

3% mi $%nj 



3 m n 



«3,,2 w í _ x 

c UJn x — ■» "rrm L'n 



mn u ij 



E 



p^m 



Sa — 3 



R(mp)< 



(A.5) 



5Z al = 3M|3w 3 E - MOW - V) 



X 



9=1 



(mq) r>(mq) r>(mq) r>{mq) ' 

-Ttj _TLj -^fc 



(A.6) 



These terms can be evaluated at equilibrium by making the following sub- 
stitutions: 



R {mn) = £S i3 {u m ~U n ), 

R^ = e\u m -u n \. 



(A.7) 
(A.8) 



We thus obtain: 
dV 



dXr 



d 2 V 



3% mi 3%nj 



MA 



•Í3 



_ sgn(u m - u n ) 
Um ^ (u - u ) 2 

n= l {""in a n) 
V n=ém 



Muj 



"mn"ij 2 / j 



2 N 



van u pn 



LO- 



l |tí m Up\ 



\ 

i 

) 



(A.9) 
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d 3 V 



3% mi dx n j dXpk 



3Mu 



2 N 

E 

9=1 
q^m 



sgn{u rn - u q )(5 mn - 5 qn )(5 mp - 8 q 



x (SijSkz + 8 ik 8 j3 + 5 jk 5 i3 - 58 a 8 j3 8 k3 ) 



(A.10) 



(A.11) 



Equation (A.9) give the equations that determine the equilibrium positions 
u m . From eq.(A.lO) we can write 



^ m,n=l í,j=l d x midx n j 

with 



and 



Mu 2 N 







m,ra=l 



* 8 -8 



p=l 

p^m 



llim. U 



l 3 : 



B 



"mn \ "ran 0p n 
^ ^p| 



(A.12) 
(A.13) 

(A.14) 



which are equivalent to the expressions given in Section 2. And, finally: 
1 A A <9 3 V 



E E 



6 m,n,p=li,j,k=l dx mi dx n jdXpk 



£,mi£,njÇp k 







E Cmnp ^p3 (2^m3^rí3 — 3^ m i^„i — 3^ m 2^n2) , (A.15) 



2£ 



mnp=l 



where 



IV 

/ii/í.p ^ ' 

q=l 
q^m 



sgn{u q - u m )(5 mn - 5 qn )(5 mp - 8 qp ) 

(u m - u q y 



(A.16) 



Appendix B 



The following tables list the non-zero vàlues of the mode cross-coupling coefi- 
cients D mnp defined by eq.(2.20) for crystals of 2 to 10 ions, together with the 
associated resonant anisotropy parameter a. Other vàlues may be determined 
via the symmetry relation D mnp = D mpn = D npm = D nmp = D pmn = D pnrn . 
Note however that the value of a is symmetric only under interchange of the 
first two indices. 
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Table 2a: Coupling constants for resonances of the first kind (in which a 
transverse phonon creates a longitudinal phonon and a transverse phonon). 



{m,n,p} 


D mn p 


Ores 


{m,n,p} 




Ores 


{m,n,p} 




Ores 


N=6 


N=7 




{3,6,3} 


2.8395e-4 


0.11575 


{4,7,3} 


5.4794e-4 


0.087700 




N=8 


{3,8,3} 


3.1152c-06 


0.062943 


{4,8,4} 


1.0031e-3 


0.069797 


{5,8,3} 


7.4469c-4 


0.068879 


N=9 


{3,9,4} 


5.7509c-06 


0.054162 


{4,9,3} 


5.7509c-06 


0.049960 


{5,9,4} 


1.2967e-3 


0.056583 


{6,9,3} 


8.8502e-4 


0.055637 














N=10 


{3,10,3} 


5.9416e-08 


0.035465 


{3,10,5} 


7.3591c-06 


0.046313 


{4,10,4} 


1.0248c-05 


0.044211 


{5,10,3} 


7.3591e-06 


0.040715 


{5,10,5} 


1.5984e-3 


0.047291 


{6,10,4} 


1.4703e-3 


0.046822 


{7,10,3} 


9.8451c-4 


0.045964 
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Table 2b: Resonances of the second kind, (in which a longitudinal phonon 
creates two transverse phonons). 



{m,n,p} 






{m,n,p} 




Ores 


{m,n,p} 






N=2 


N=3 


{2,2,2} 


-1.1225 


0.57143 


{3,2,3} 


-1.5754 


0.30917 


{3,3,2} 


-1.5754 


0.31746 


N=4 


{3,3,4} 


-1.8332 


0.21132 


{4,2,4} 


-1.9493 


0.19337 


{4,3,3} 


-1.8332 


0.20560 


{4,4,2} 


-1.9493 


0.20391 


{4,4,4} 


2.1084 


0.15429 








N=5 


{4,3,5} 


-1.9683 


0.14895 


{4,4,4} 


-0.89355 


0.15387 


{5,2,5} 


-2.2887 


0.13340 


{5,3,4} 


-1.9683 


0.14187 


{5,4,3} 


-1.9683 


0.14619 


{5,4,5} 


3.1611 


0.11561 


{5,5,2} 


-2.2887 


0.14311 


{5,5,4} 


3.1611 


0.11668 








N=6 


{4,4,6} 


-1.8850 


0.11439 


{5,3,6} 


-2.0625 


0.10985 


{5,4,5} 


8.1287e-2 


0.11527 


{5,5,6} 


4.2528 


0.092393 


{6,2,6} 


-2.6084 


0.098229 


{6,3,5} 


-2.0625 


0.10398 


{6,4,4} 


-1.8850 


0.10784 


{6,4,6} 


4.0185 


0.088947 


{6,5,3} 


-2.0625 


0.10937 


{6,5,5} 


4.2528 


0.091510 


{6,6,2} 


-2.6084 


0.10658 


{6,6,4} 


4.0185 


0.091151 


{6,6,6} 


-2.3706 


0.075762 














N=7 


{6,3,7} 


-2.1382 


0.084291 


{6,5,7} 


4.9336 


0.074398 


{6,6,6} 


4.3077 


0.075526 


{7,2,7} 


-2.9149 


0.075728 


{7,3,4} 


5.4794e-4 


0.088234 


{7,3,6} 


-2.1382 


0.079773 


{7,4,5} 


-1.7957 


0.082749 


{7,4,7} 


4.7748 


0.070376 


{7,5,4} 


-1.7957 


0.084569 


{7,5,6} 


4.9336 


0.072836 


{7,6,3} 


-2.1382 


0.085059 


{7,6,5} 


4.9336 


0.074000 


{7,6,7} 


-3.8751 


0.062565 


{7,7,2} 


-2.9149 


0.082796 


{7,7,4} 


4.7748 


0.073152 


{7,7,6} 


-3.8751 


0.062861 














N=8 


{6,6,8} 


5.2871 


0.062147 


{7,3,8} 


-2.2035 


0.066792 


{7,5,8} 


5.4173 


0.060792 


{7,6,7} 


3.6190 


0.062378 


{7,7,6} 


3.6190 


0.062674 


{7,7,8} 


-5.8640 


0.053277 


{8,2,8} 


-3.2119 


0.060393 


{8,3,5} 


7.4469e-4 


0.069615 


{8,3,7} 


-2.2035 


0.063341 


{8,4,6} 


-1.7158 


0.065623 


{8,4,8} 


5.4691 


0.057068 


{8,5,5} 


-1.5767 


0.067200 


{8,5,7} 


5.4173 


0.059134 


{8,6,4} 


-1.7158 


0.068081 


{8,6,6} 


5.2871 


0.060526 


{8,6,8} 


-5.1946 


0.052159 


{8,7,3} 


-2.2035 


0.068170 


{8,7,5} 


5.4173 


0.061000 


{8,7,7} 


-5.8640 


0.053018 


{8,8,2} 


-3.2119 


0.066389 


{8,8,4} 


5.4691 


0.060029 


{8,8,6} 


-5.1946 


0.052896 


{8,8,8} 


2.1526 


0.046043 








N=9 


{7,6,9} 


5.4146 


0.052246 


{7,7,8} 


1.6006 


0.053125 


{8,3,9} 


-2.2623 


0.054311 


{8,5,9} 


5.7898 


0.050462 


{8,6,8} 


2.6715 


0.052000 


{8,7,7} 


1.6006 


0.052867 


{8,7,9} 


-7.3423 


0.045500 


{8,8,6} 


2.6715 


0.052744 


{8,8,8} 


-7.5531 


0.045900 


{9,2,9} 


-3.5017 


0.049433 


{9,3,6} 


8.8502c-4 


0.056425 


{9,3,8} 


-2.2623 


0.051650 


{9,4,5} 


1.2967c-3 


0.056883 


{9,4,7} 


-1.6468 


0.053425 


{9,4,9} 


6.1214 


0.047245 


{9,5,6} 


-1.4030 


0.054734 


{9,5,8} 


5.7898 


0.048927 


{9,6,5} 


-1.4030 


0.055586 


{9,6,9} 


-6.3927 


0.044000 


{9,7,4} 


-1.6468 


0.056011 


{9,7,6} 


5.4146 


0.050996 


{9,7,8} 


-7.3423 


0.044970 


{9,8,3} 


-2.2623 


0.055954 


{9,8,5} 


5.7898 


0.051126 


{9,8,7} 


-7.3423 


0.045400 


{9,8,9} 


3.8083 


0.039874 


{9,9,2} 


-3.5017 


0.054558 


{9,9,4} 


6.1214 


0.050180 


{9,9,6} 


-6.3927 


0.045081 


{9,9,8} 


3.8083 


0.039988 


N=10 


{7,7,10} 


5.2107 


0.044976 


{8,6,10} 


5.4385 


0.044332 


{8,7,9} 


-0.42626 


0.045380 


{8,8,8} 


-1.9896 


0.045779 


{8,8,10} 


-8.6478 


0.039673 


{9,3,10} 


-2.3165 


0.045099 


{9,4,9} 


3.3150 


0.047113 


{9,5,10} 


6.0930 


0.042512 


{9,6,9} 


1.6336 


0.043878 


{9,7,8} 


-0.42626 


0.044846 


{9,7,10} 


-8.5012 


0.039123 


{9,8,7} 


-0.42626 


0.045279 


{9,8,9} 


-8.0166 


0.039760 


{9,9,6} 


1.6336 


0.044957 


{9,9,8} 


-8.0166 


0.039874 


{9,9,10} 


6.3437 


0.035161 


{10,2,10} 


-3.7855 


0.041300 


{10,3,7} 


9.8451 e-04 


0.046733 


{10,3,9} 


-2.3165 


0.043017 


{10,4,6} 


1.4703c-3 


0.047276 


{10,4,8} 


-1.5874 


0.044422 


{10,4,10} 


6.7435 


0.039798 


{10,5,7} 


-1.2642 


0.045500 


{10,5,7} 


-1.2642 


0.045500 


{10,5,9} 


6.0930 


0.041163 


{10,6,6} 


-1.1721 


0.046259 


{10,6,8} 


5.4385 


0.042265 


{10,6,10} 


-7.5034 


0.037564 


{10,7,5} 


-1.2642 


0.046721 


{10,7,7} 


5.2107 


0.043067 


{10,7,9} 


-8.5012 


0.038483 


{10,8,4} 


-1.5874 


0.046923 


{10,8,6} 


5.4385 


0.043500 


{10,8,8} 


-8.6478 


0.039085 


{10,8,10} 


5.3952 


0.034687 


{10,9,3} 


-2.3165 


0.046823 


{10,9,5} 


6.0930 


0.043468 


{10,9,7} 


-8.5012 


0.039269 


{10,9,9} 


6.3437 


0.035057 


{10,10,2} 


-3.7855 


0.045724 


{10,10,4} 


6.7435 


0.042604 


{10,10,6} 


-7.5034 


0.038860 


{10,10,8} 


5.3952 


0.035000 


{10,10,10} 


-1.7362 


0.031318 
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